Apparent Hubble acceleration from 
large-scale electroweak domain structure 
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The observed luminosity deficit of Type la supernovae (SNe la) at liigfi redsliift z can be explained 
by partial conversion to weak vector bosons of photons crossing large-scale electroweak domain 
boundaries, making Hubble acceleration only apparent and eliminating the need for a cosmological 
constant A > 0. 
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I. INTRODUCTION 

After the initial surprise caused by the announcement 
of a luminosity deficit from high redshift z Type la super- 
novae (SNe la) dii, the consensus quickly emerged 
that we are witnessing an accelerating Hubble expansion, 
implying a cosmological constant A > 0. This has forced 
a profound change in our view of the large scale struc- 
ture and composition of the universe, leading us from a 
preferred model with negative curvature and density pa- 
rameter n « 0.3 dominated by (mostly dark) matter to 
a flat geometry with 17 « 1 dominated by the dark en- 
ergy term IIa ~ 0.7 (ACDM). Most importantly from the 
point of view of fundamental physics, it has also left us 
with the massive embarrassment of a finite A 120 orders 
of magnitude below the Planck scale, where it's custom- 
arily argued on dimensional grounds that the zero point 
energy of quantum fields coupled to gravity should show 
up, absent some fundamental principle forcing it to van- 
ish exactly. 

These far-reaching implications have motivated many 
studies of alternative explanations for the dimming of 
high-z SNe la, from variations in intrinsic luminosity 
(chemical abundances, stellar populations) through mod- 
ified light propagation (gravitational lensing, gray dust) 
to systematic observational errors (selection bias) 
But to date, all suggested mechanisms have proved inca- 
pable of producing effects of the needed size, -0.25 mag- 
nitudes at z ~ 0.5 (where the luminosity deficit has its 
greatest leverage on A), i.e. a ratio between observed 
luminosity I and expected luminosity Ie 
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The case for dark energy rests squarely on this number. 
In spite of common claims to the contrary, dimming of 
SNe la (and now of gamma ray bursts [7|) remains the 
only direct evidence of accelerated expansion Q. Even 
disregarding the possibility of significant systematic er- 
rors in distance measures the oft-quoted consis- 
tency of ACDM with other observations is only a neces- 
sary condition for its validity, not a sufficient one. The 



true power of this condition is its ability to falsify the 
framework: find an explanation for one data set leav- 
ing insufficient room for accelerated expansion to fit the 
others and the entire framework is invalidated. 

The main objective of this paper is to show how the 
standard electroweak and big bang models can lead to the 
observed ^|lE without any Hubble acceleration actually 
taking place, and how they are in fact constrained by 
observation to do so in a way which leaves little room for 
A > 0. 



II. PHOTON TRANSFORMATIONS 

In the standard SU(2)xU(l) model of electroweak in- 
teractions [iH[i2|[l3' the U(l) gauge field B^{x) and the 
three SU(2) gauge fields W^^(x) are collected in the ma- 
trix (the field-dependent part of the covariant derivative) 



g'B^ {x) + gWl (x) gW^ (x) - igW^ [x) 
gWl,{x) + igWlix) g'B^{x) - gW^x) 



(2) 



[g — SU(2) couphng constant; g' = U(l) coupling con- 
stant) acting on weak isospin doublets. Given an arbi- 
trary 7W^i(x), we can decompose it into individual gauge 
fields using 
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where all traces are understood to be over weak isospin 
space only, tq = I2 and ti, T2, T3 are the Pauli matrices 



(7) 



The photon A^j^(x) and the neutral weak vector boson 
Z^{x) are defined as the linear combinations 
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cos(6'w ) 
sin(6'vt') 
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where 9w is the Weinberg (weak mixing) angle, 

gsa\{6w) = g' cos{6w) = e > (9) 

(e = electric charge of the proton). Combining ([6]) 
and ([8]), we can therefore write 



sm[Uw ) 

cos{6w) 
25' 



Ti-[t3M^{x)] 
Tr[roM^(a;)] 
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Setting Z°{x) = W^{x) = W^l^x) = and A^,{x) ^ 
1, inverting ([8]) and plugging the results into ^ yields 
M^{x) for a normalized pure photon state: 
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(completely delocalized and therefore monochromatic; an 
envelope can be imposed without consequence for the 
present argument). Now consider the effect on M^{x) of 
a global SU(2)xU(l) transformation 



Z 2 
U = CXp ( - CJoTo + - LUjTj 



(12) 



with U(l) parameter and SU(2) parameters u) = 
(tJi,LJ2,^3): 

M^{x) ^ Ml{x) = UM^{x)U^ (13) 

Substituting this into (jlOp and using the cyclic property 
of the trace, 



--(^>^).Tr[[/tr3C/M,(x)] 

■Tr[M^(x)] (14) 
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U(l) transformations associated with woTq drop out. 
Specializing to the pure photon state of (jlip . in- 
troducing the shorthand 

cj^ = wi+w| (15) 
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and doing the traces, (fT4|) reduces to 



with 



(16) 
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(a;) = sin^' 



cos(u;) + wf 
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(see Fig. [H O. Together with (dTj), equation (dH) de- 
scribes the effect of a global transformation with SU(2) 
parameters u) (and arbitrary U(l) parameter coq) on a 
pure photon state. On its own, since A^{x) is propor- 
tional to the photon number operator, it gives us the 
fraction of photons surviving the transformation (the rest 
having turned into Z'^s and linear combinations of W^s 
and W^^s, i.e. W^s). The full import of this residual 
luminosity £{u)) will become evident in Section IVlII 
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FIG. 1; Residual luminosity £{(jj) for ij± G [0, \/327r], cus £ 
[0, 47r] (contour plot). 




FIG. 2; Residual luminosity £{lv) for uj± G [0, V327r], ws G 
[0, 47r] (perspective plot). 



III. ELECTROWEAK PHASE TRANSITION 



In order to give the standard model fields (inertial) 
mass without explicitly breaking the gauge symmetry 
of the Lagrangian (needed for renormalizability [l3|[lB|)j 
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the weak isospin doublet 



is introduced, where 4)^{x) and (jP^x) are complex 
Lorentz scalars transforming under SU(2)xU(l) accord- 
ing to 

= U<^{x) (20) 

with the U of (fT2)) . Augmenting the Lagrangian with a 
potential term featuring a manifold of degenerate ground 
states for <^( x) ^ (all connected by U) turns ^{x) into 
a Higgs field (l7| which breaks the symmetry dynam- 
ically by picking a non-zero VEV (vacuum expectation 
value) ($) in a specific ground state. The asymmetric 
VEV splits the interaction terms between $(a;) and the 
other fields into effective mass terms and new interac- 
tion terms involving <^(x) excitations above the selected 
ground state^. 

In the early universe, this electroweak phase transition 
(EWPT) is expected to have occurred when the temper- 
ature of the expanding plasma fell below the electroweak 
scale. Tew ~ 250 Ge^ « Plugging T^vf into 

the standard big bang model yields a cosmological scale 
factor ~ 10~^^ and a Hubble radius rn ~ 2.9-10"^^ rela- 
tive to the present epoch, i.e. causally connected patches 
roughly 1 cm across, each selecting a VEV independently 
of the others. 

Much effort has gone into studying the onset and early 
stages of the EWPT, primarily due to their implications 
for baryogenesis fToj [20| . Scenarios motivated by the lat- 
ter generally begin with 0(1) bubbles of ($) ^ nucleat- 
ing within each causal patch and then expanding into the 
surrounding plasma at some small fraction of the speed 
of light (typically ^ 0.1; but see [2l[ for much lower 
estimates) determined by the equilibrium between out- 
ward pressure and plasma friction. Colliding bubbles are 
believed to have undergone repeated bounces and reheat- 
ing before finally merging [22| , leading to a period of slow 
average growth lasting several orders of magnitude longer 
than their initial expansion. Such a (first order, i.e. bub- 
bly) phase transition has been ruled out for the minimal 
standard model 23] [24] but remains a viable possibility 
in its various extensions. 

In the minimal standard model, the EWPT is a con- 
tinuous crossover, i.e. an intrinsically non-perturbative - 
and therefore analytically challenging - process. While 
one may argue that the natural speed scale in this case 
is that of sound in a relativistic plasma, Vg — ^1/3, lit- 
tle is actually known with certainty about the dynamics 
involved. 



For an alternative, perhaps more intuitive derivation of iijjS), 
diagonalize the gauge boson mass matrix for arbitrary (<&) and 
take the scalar product of massless eigenstates, i.e. photons, for 
VEVs related by the global transformation parameter a> [iSll . 



IV. VACUUM REALIGNMENT 

Fortunately, apart from an overall scale factor set by 
the amount of time needed to convert all space to the bro- 
ken symmetry phase, the late epoch should be relatively 
unaffected by the early dynamics. Once all space has 
been converted, the problem boils down to the realign- 
ment of ($) across adjacent causal patches, no matter 
how they got their initial VEVs. 

To avoid a common source of confusion, we need to 
clarify the meaning of the term "realignment" . In the 
absence of gauge interactions, it is unambiguous: when 
'^{x) varies across space, there is gradient energy, which 
is minimized by $(2:) evolving to a constant (<&). With 
gauge interactions, as long as ^{x) does not go to zero 
anywhere in the volume under consideration, it is al- 
ways possible to substitute a position-dependent U{x) 
into ((20|) and "gauge away" any misalignment in ^{x). 
But such a local transformation requires that we also use 

M^{x)^M'^{x) = -iU{x)[d^U\x)\ 

+ U{x)M^,{x)U\x) (21) 

in lieu of What this amounts to is a change of 

variables: we are trading gradients in (f>(x) for gradients 
in M^,{x) 2. 

The obvious way to cure the resulting ambiguity is to 
choose a (complete) gauge fixing condition and then stick 
with it. The natural choice when considering multiple 
vacua (as opposed to doing perturbation theory in one 
of them, the traditional business of particle physics) is 
a physical gauge which leaves the internal symmetry of 
the theory manifest, allowing the concept of alignment 
to retain its intuitive meaning (unless otherwise stated, 
this point of view will be implied for the rest of this 
paper). But the ultimate arbiter is always energy. A 
true vacuum is a global energy minimum of the theory; if 
the energy density within a given volume is everywhere at 
such a minimum, we say that we have alignment within 
that volume, even though our particular choice of gauge 
may make individual field components look anything but 
aligned. 

This brings us to another, closely related source of con- 
fusion: the gauge trajectory argument. If moving along a 
spacetime trajectory takes us through a sequence of $(a;) 
values related by a gauge trajectory, i.e. by a sequence of 



^ Introductory QFT texts tend to focus on infinitesimal transfor- 
mations, sometimes leaving the full form II21I I to more advanced 
discussions of non-perturbative solutions. Unfortunately, this 
practice seems to have spawned a legion of phenomenologists who 
genuinely believe that the standard electroweak model has only 
one physical vacuum, to which all other vacua can be transformed 
without ulterior consequences. What they make of sphalerons, 
the saddle point solutions connecting different vacua which un- 
derpin the relevance of electroweak theory to baryogenesis, is a 
mystery. 
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successive infinitesimal transformations, then (and only 
then) using the inverse of this gauge trajectory for U (x) 
in ((20|) and ([2T|) will realign <I>(a;) along the spacetime tra- 
jectory without affecting the energy carried by the gauge 
fields. The oft-quoted argument, due to Turok, that tex- 
tures (knot-like, gradient-only <i>(x) configurations) be- 
come true vacuum configurations when $(2:) has gauge 
interactions [1^ is a special case of this observation. It 
does not imply, as it's sometimes misconstrued to do, 
that gradient-only <I>(a;) configurations residing entirely 
on the vacuum manifold are trivial in gauge theories, 
only that they can generally be expected to be unsta- 
ble (Turok's use of the word "become" may be partly to 
blame for this common misunderstanding; it should be 
read as "evolve to"). If we arrange $(a:) in such a con- 
figuration on a spacelike 3-surface with all gauge fields 
M^(x) set to zero, the result is indistinguishable from 
having the same configuration in the corresponding un- 
gauged theory, and so will necessarily have positive en- 
ergy. If we now let the field equations run their course, 
we will see the energy being dispersed as 7V/^(a;) picks 
up. There is no question about this being a real, physi- 
cal process playing out over time. The only question is: 
how much time?'^ 

In the case of a localized $(a;) configuration living in 
an otherwise empty vacuum (such as a texture), energy 
can disperse in all directions at the speed of light; in 
practice the only relevant time scale is that of the gauge 
interactions. The natural expectation is then for the con- 
figuration's peak energy density to decay exponentially 
with a half-life on the interaction scale (see [23] for actual 
simulation results on electroweak texture decay). In the 
more complicated case of a random configuration filling 
all space, there is a second time scale: the average prop- 
agation time to the first recurrence of the field values 
(derivatives included), better known as inverse tempera- 
ture. Once we hit such a recurrence - and in an infinite 
space, one is guaranteed to occur to any desired degree of 
precision in every spatial direction - we have a periodic 
boundary condition, implying conservation of energy, as 
opposed to the absorbing boundary conditions which al- 
low localized configurations living in an empty vacuum 
to simply vanish from sight. The natural expectation 
is then for energy density to settle into a semi-periodic 
pattern. 

The weak spot in the gauge trajectory argument should 
now be evident: the argument tells us that M'^(x) can be 
arranged so that a ^(x) configuration which never leaves 
the vacuum manifold may be gauged away "for free" , 
but it does not relieve us of having to explain where the 
energy of the initial field configuration ends up, nor how 
it is carried away. No explicit mechanism, no decay. 



^ Turok argued that the answer is a microphysical time, set by 
the time scale of gauge interactions, but the argument fails when 
massless gauge fields remain after SSB. See |26i | and Sections 
IV. A, IV.B in O. 



It is therefore necessary to consider the flow of con- 
served quantities between ^{x) and all other fields, 
fermions included. In this paper, we are primarily con- 
cerned with $(a;) configurations which can interact with 
photons, i.e. which carry electric charge. A single, mas- 
sive, charged electroweak boson - the ultimate localized 
state of the theory - will of course decay on the time 
scale of electroweak interactions, but that's beside the 
point; the issue at hand is the lifetime of extended field 
configurations on the vacuum manifold, not of isolated 
excitations above it. To state the obvious, a boson VEV 
is not a "dust cloud" of distinct on-shell particles (if it 
were, the electroweak vacuum itself, being a Higgs con- 
densate, would decay as fast as a lone Higgs boson); it's 
a continuous quantity whose evolution is determined by 
the field equations derived from the theory's effective ac- 
tion (classical action + radiative corrections). 

Fortunately, we need not carry out the full program 
here. Thanks to the work on pai r production initiated 
a long time ago by Schwinger [28|, we know that it can 
be viewed as a tunneling process putting virtual fermion 
pairs on shell at the boson field's expense. As long as 
a spatial volume subtended by the Compton wavelength 
l/rUe of the lightest charged fermion, the electron, con- 
tains at least unit electric charge and gradient energy 
> TTie « 511 keV (plus the negligible rest mass of a neu- 
trino) there can be spontaneous emission of electron -|- 
antineutrino pairs. When charge and gradient energy 
densities fall below these thresholds, the tunneling rate 
becomes exponentially suppressed. (The elementary in- 
sight that charge conservation limits the decay rate of a 
charged boson condensate, making it quite distinct from 
the decay rate of a single boson, can also be arrived at 
by traditional statistical mechanics, as demonstrated to 
dramatic effect in 29]). Further dissipation must then 
be catalyzed by interactions with the environment. 

In cosmology, interactions are strictly between field 
modes with wavelength < ru- Short of renouncing local- 
ity, superhorizon modes are decoupled. No local interac- 
tion can therefore dissipate conserved quantities carried 
by such modes. If interactions between a field and the 
environment effectively shut down at some point in time, 
field modes which were outside the horizon at that time 
are still around today. 

Consider a ^{x) configuration looking to lose some pos- 
itive electric charge. The cheapest catalyst is a free elec- 
tron. Turning it into a 511 keV neutrino does not cost 
any energy, but a spatial charge density > e is still 
required (the matrix element is the same as for pair cre- 
ation, we have simply reversed an external momentum). 
Once charge density falls below this threshold, it can no 
longer be dissipated away effectively, and so neither can 
the ^(x) configuration carrying it. Since the (l)'^{x) has 
unit electric charge, this implies a threshold energy den- 
sity ~ mHrnl = m'^{me/rnH)^, corresponding to a tem- 
perature factor ~ {nie/mHY^^ ~ 10^^ relative to Tew^ 
i.e. T ~ lOMeV « lO" K; the leptonic era. Relative to 
the present epoch, the cosmological scale factor was then 
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~ 10~^° and the Hubble radius 10~^^, or 10^ km. 
This gives us the scale beyond which $(2:) configurations 
should have been safe from dissipation. 

Incidentally, the subsequent redshifting of 10 AleV by 
the scale factor 10~^° lands us right at today's dark en- 
ergy scale, 10^"^ eV, suggesting that Higgs gradients may 
provide the missing energy density required for SI « 1. 

To state the obvious once more, 10"^ km is a macro- 
scopic distance separated by some 23 orders of magnitude 
from the electroweak scale, the natural focus of works on 
the EWPT. There is therefore no contradiction between 
the commonly expected fast dissipation of configurations 
with characteristic sizes on the electroweak scale, like Z 
strings (unless stabilized by plasma effects, as argued by 
Nagasawa and Brandenberger [s^l) and long-lived modes 
with wavelength > 10'^ km. 

Summing up, the misalignment in ($) after the EWPT 
is a physical reality which can not simply be "gauged 
away" . As the Hubble radius grows, regions establishing 
causal contact for the first time after the transition have 
to realign according to the field equations. This realign- 
ment can not proceed at a speed faster than light's. 

A final clarification has proved necessary: 

In theories supporting nontrivial mappings between 
spacetime and internal symmetry space, realignment 
sooner or later hits a stopping point in the form of a 
topological defect, i.e. a (<i>) configuration which can 
not be realigned throughout all space without leaving 
the vacuum manifold, at an energy cost on the order of 
the symmetry breaking scale (Kibble mechanism [3l|). 
Below this scale, such defects are therefore classically 
stable. Their importance for cosmology has been well 
understood for decades and has spawned a vast litera- 
ture. In particular, it has long been appreciated that 
domain walls, i.e. topologically stable, two-dimensional 
defects passing through ($) — 0, would be catastrophic, 
as each Hubble volume would quickly become dominated 
by a massive, single wall [s^llljl. Theories with domain 
wall solutions are therefore ruled out by observation. 

The domain structure referred to in the title 
has nothing to do with topologically stable do- 
main walls. The term "domain boundary", as opposed 
to "domain wall" , was adopted to help keep this distinc- 
tion in mind. 



V. DOMAIN STRUCTURE 

Unlike the onset and early stages of the EWPT, its 
later stages have attracted little attention. Since the 
standard electroweak model features no topological de- 
fects and no (known) dynamically stable solutions ^, to 
the extent that realignment has even been recognized as a 



* A stable oscillating solution (and potential non-exotic dark mat- 
ter candidate) was found numerically after this was written |34| . 



real, physical process, it has simply been assumed to pro- 
ceed at the rate typical of electroweak interactions, with- 
out any consequences worthy of notice. There is therefore 
little in the way of past results to help us gain some in- 
sight into its dynamics. Ultimately, settling the issue will 
come down to massive numerical simulations. Until such 
studies are performed, the analytical intractability of the 
highly non-linear electroweak field equations forces us to 
rely on heuristics and on analogy with other physical sys- 
tems. 

In Section IIVI conservation of energy led us to expect 
the emergence of a semi-periodic spatial pattern. A fur- 
ther inroad to the problem is again provided by conser- 
vation of electric charge (Q). 

Fixing a global coordinate system in weak isospin space 
implies fixing a definition for Q. By the standard con- 
ventions, the two complex Higgs components 4''^{x) and 
4P{x) of carry unit and zero Q, respectively. While 
the universe as a whole is assumed to be electrically neu- 
tral, the random choice of ($) at the EWPT will therefore 
initially result in a random charge distribution. 

Consider a volume with radius R ^ rn- It starts off 
containing some net Q. On average, this net charge will 
flow outward. Since the initial distribution is random, 
the trajectory of a charged test particle starting from 
the center of the volume is a three-dimensional random 
walk, covering an average distance 

(r) cx Vi (22) 

in time t and becoming fractal-like (i.e. statistically scale 
invariant) at late times. When (r) = R, the distribu- 
tion has become a set of charged shells enclosing neutral 
regions with average radius R. As they collide, oppo- 
sitely charged shells annihilate, their enclosed volumes 
merging; equally charged shells press against each other. 
By comparison with compression of randomly packed 
spheres, the resulting partitioning of space can be ex- 
pected to consist of irregular polyhedra with an average 
of w 13 faces [Hd^. 

The upshot is that the current operator for a complex 
scalar field is all derivative terms. By tracking Q flows, 
we are therefore tracing out regions where ($) is not con- 
stant. These are our domain boundaries. As time goes 
by and domains merge, average domain size grows, re- 
ducing the number of domain boundaries and their total 
surface area. 

What we have here are all the essential features of a 
problem well known to materials scientists: local inter- 
actions causing the formation of flat surfaces separating 
polyhedral domains, followed by minimization of total 
surface area by successive domain merging. Such a sys- 
tem is known as a foam, the problem of its evolution as 
"foam coarsening" (or "grain growth"). Its main attrac- 
tion is universality: the underlying interactions do not 
matter as long as they provide the above features. Poly- 
hedral solutions are indeed commonplace in non-linear 
fleld theories and have also been explicitly constructed 
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FIG. 3: Two semi-regular tetrakaidecahedra superimposed 
on an illustration of their section from Kelvin's paper On 
the Division of Space with Minimum Partitional Area (1887). 
Real foams do not display such regularity, but their average 
properties resemble those of the tetrakaidecahedron. 



in the standard electroweak model^ WA - 

Foam coarsening is a tricky problem, von Neumann fa- 
mously solved the two-dimensional case on the spot upon 
hearing about it, equating the area rate of change of a 
two-dimensional domain to the number of its sides [38| . 
but the three-dimensional case remains an open chal- 
lenge. Two universal rules, first described by Plateau in 
his classic work on soap films (39| . are known to apply: 
along each edge, three faces of the constituent polyhedra 
meet at angles of 27r/3; at each vertex, four edges meet 
at the tetrahedral angle arccos(— 1/3). As a consequence, 
the average number of faces (/) and the average number 
of edges per face (rie) can be shown to satisfy 
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6 - (ne 



(23) 



The fundamental difficulty in going beyond this level 
- and the crucial point of this section - is that while the 
shaping of domain boundaries is a local process proceed- 
ing on the time scale of the underlying interaction, the 
foam's evolution, i.e. domain growth, is not. When two 
domains merge, it is not an isolated event: their nearest 
neighbors (and then their nearest neighbors, and so on) 
are also affected, ultimately forcing the whole foam to 
rebalance. In the case at hand, currents will speed up, 
slow down or even reverse, taking our charged test parti- 
cle along on a continued three-dimensional random walk. 
Since currents live on domain boundaries, we are led to 



^ Better still, the low energy effective field theory of the elec- 
troweak boson sector is a gauged non-linear sigma model 
(NLSM), which is easily verified to admit the large number of 
solutions known from the plain NLSM, including polyhedral ones 



expect the statistical properties of the test particle's tra- 
jectory, i.e. scale invariance and average radius growth 
according to (f^^ . to carry over to the foam itself. 

This does indeed turn out to be the case. After much 
experimental and numerical work (not least on Potts 
models, believed to closely approximate the finite tem- 
perature dynamics of SU(N) gauge theories in general 
and of SU(2) in particular [40|), a consensus has emerged 
in recent years that the average volume growth rate of 
three-dimensional domains with / faces is well described 
by a simple linear dependence on / [4l| [i^l 



{Vf) 



-1/3 



d{Vf) 
dt 



= «(/-/o) 



whence 



/o 



if) 



(24) 



(25) 



Empirically, /q ~ 14±2, well in line with our guesstimate 
based on sphere crunching. Growth laws on the form (j24p 
(i.e. with right hand side depending only on the number 
of faces, or more generally on topological features of the 
domains) lead to the anticipated time dependence ([^^ 
for the average domain radius. This result is also ex- 
pected on dimensional grounds from the emergence of a 
scaling state, characterized by growing average domain 
size but time-independent topological and area distribu- 
tions [43] , again consistent with expectations from the 
random walk argument (scale invariance). 

The scaling state is a disordered one, not a regular 
structure composed of individually near-optimal parti- 
tions such as Kelvin's famous 14-faced tetrakaidecahe- 
dron (Fig. [3|). The discovery by Weaire and Phelan 
[Hi that an 8-cell "repeat unit" with (/) = 13.5 and 
(rie) = 5.11 achieves a more efficient (smaller total sur- 
face) partitioning of a given volume than Kelvin's solu- 
tion explains why: a combination of many cells of dif- 
ferent shape can actually have lower total surface energy 
than a regular structure. 

It bears emphasizing that these are universal results, 
depending only on the assumption that the system strives 
for surface area minimization, not on any details of the 
underlying interactions. In equation p4|) . all interaction 
dependence is encapsulated in the constant k. In partic- 
ular, it can not be underscored enough that the domain 
boundaries referred to here have nothing in common with 
topologically stable defects: they do not depend on non- 
trivial mappings between spacetime and internal degrees 
of freedom for their existence, they have low energy den- 
sity, they give rise to a completely different domain struc- 
ture (the foam) , and they positively must be unstable for 
domain growth and emergence of the scaling state to be 
at all possible. Domains grow by merging, i.e. by the 
disappearance of domain boundaries; if the latter were 
stable, growth could not happen. The separation of time 
scales and the decelerating evolution embodied in ([^^ 
are collective (and essentially geometric) properties of the 
whole foam, not of individual domain boundaries. 
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VI. NUMBER OF DOMAINS 



In a relativistic setting, ([221) must be qualified by the 
requirement that the speed of Ught not be exceeded (lo- 
cally). The simplest ansatz satisfying this condition is 



(26) 



valid for t >— ta, where ta marks the onset of the scaling 
state, 



{r{ta)) = a 



with initial growth rate 



and asymptotic growth rate 



lim {f{t)) = ^ 



(27) 



(28) 



(29) 



Putting ([26|) on the standard Robertson- Walker metric 
with scale factor a(t) yields the proper ensemble domain 
radius 



{rRw{t)) 



aft) 
a(ta) 



it)) 



dr 



a{ta) 



with 



m) 



a/32 



2^1 + /32(t-i„) 



(r(r)) (30) 



(31) 



Sweeping (3 over the range [0, 2/a] will take us through 
the possible values of , the average number of observ- 
able ($) domains in an arbitrary spatial direction. What 
can we expect to find? 

The limit case /3 = is easily understood: it describes 
domains of average size a being effectively "frozen out" 
at t — ta (by their growth rate becoming negligible com- 
pared to that of the universe) and then simply coasting 
along with their comoving volume. Assuming flatness, 
if a is the Hubble radius at the time of the EWPT, the 
opposed effects of subsequent horizon growth and metric 
expansion work out to TV* ~ 2^q23-15 _ ^qS^ freeze- 
out occurs later, N,^ can be substantially lower. Section 
IIVI suggests iV$ < 10^. One frozen domain per Hubble 
radius at recombination, z « 1089 according to cosmic 
microwave background (CMB) data, would translate to 
only ~ 30 today. 

At the other end of the /3 range, the Hubble radius 
poses an absolute limit. As recently emphasized by Pen- 
rose [i^ in a revival of the old homogeneity problem of 
pre-inflationary cosmology, light reaching us now from 
quasars in opposite directions must have originated in 
different electroweak domains, as those sources have not 




FIG. 4: Expected domain boundary redshifts for a/rnita) G 
[1, 10], {r{ta)) = 10"*, ta = t{za = 1089) (recoHibination) in 
a flat, dust-matter Friedmann model. 



been in causal contact since before the EWPT, which 
occurred well after the end of inflation ^. Thus, 



1 < < 10^ 



(32) 



As will become evident in Section fVIH too small a num- 
ber, i.e. a rate of realignment too close to the speed of 
light, would indeed be incompatible with the observed 
isotropy of the universe. 

For a more detailed view, we must turn to equation 
([50)1 . By isotropy, we should expect an infinite sequence 
of domain boundaries at (average) proper distances 



dn{t)^{rRw{t)){2n-l) ng[l,( 



(33) 



What we actually observe is light emitted at time te 
and reaching us at time t, so 



dr 



- dji (tern ) 



(34) 



Given a{t), i.e. a cosmological model, we can solve ([M)) 
for tem and plug the result into 



aft) 



1 



(35) 



^ A word of caution is in place here. As it stands, Penrose's ar- 
gument [45I I can be interpreted as building on the notion that 
the Weinberg angle 9^ was chosen randomly at the EWPT (see 
pages 651 and 743), just like the direction of (<I>). While this 
may be the case in extended theories, it is not how the standard 
electroweak model works. In the standard model, S^v j^^t 
a parameter, not a field, and its value is the same throughout 
spacetime. 
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space 



(36) 



FIG. 5: Probability distribution functions Pn{(-) (vertically 
rescaled) for the residual luminosity £ after n boundary cross- 
ings between n + 1 domains with independent ($) . 



to obtain the observed redshift z. In practice, this must 
be done numerically. See Appendix [A] for easily adapt- 
able R [i^ code and Fig. Ufor a simple example (but not 
necessarily an unrealistic one, at least for small z): a flat 
dust-matter model featuring 0(1) slow-growing domains 
per Hubble radius at recombination. 

To recap the story so far: while average domain radius 
is cx the Hubble radius is cx t (assuming flatness). 
Therefore, while any finite volume will eventually find 
itself within a single domain, the number of observable 
domains grows with time (cx V7, again assuming flat- 
ness) . 



VII. RESIDUAL LUMINOSITY 
DISTRIBUTIONS 



Imagine a physicist in the (<&) domain !Di preparing a 
pure photon state as defined by (fTTj) and sending 
it to a colleague in the adjacent domain 2)2. To verify 
the purity of the received 7W^, the colleague performs 
the operation (jlOp on it (e.g. by measuring the peak 
interference amplitude with a locally produced reference 
state). What will she see? 

In their work, both physicists must implicitly rely on 
their respective, local definition of a photon. The obvious 
way to distinguish a photon from a is that only the 
latter has mass, a fact which is made explicit in the elec- 
troweak Lagrangian by a change of basis in weak isospin 



making € 3? the only non-zero component of ($). This 
choice of basis aligns our coordinate axes in weak isospin 
space with the direction which was picked by the Higgs 
field at the EWPT. The SU(2)xU(l) (sub)symmetry of 
the ground state manifold guarantees the existence of a 
transformation U, as defined by ([T^ . which achieves this 
alignment. The Higgs and gauge fields in the new basis 
are given by ([20|) and ^3]). respectively. In this basis, 
the has a mass term and the photon does not. 

But a global SU(2)xU(l) transformation U can put 
($) on the form ([36]) in only one of the two domains Di 
and 2)2. In the other domain, where ($) is different, 
U has no particular significance. Therefore, the local 
definition of a photon is not the same in Si and in ©2^- 

Formally, it is of course possible to perform a local 
transformation and "gauge away" any misalignment in 
($) everywhere in the universe, but as we saw in Section 
rVl only at the cost of using (pij) in lieu of i.e. of 



trading gradients in ($) for gradients in Mf^{x). This 
amounts to a change of variables and of focus, from the 
Higgs field inside domains to electrically charged gauge 
fields across their boundaries. While it may facilitate the 
description of processes at the boundaries, it does not 
affect the physics, and it does not help our two physicists 
somewhere inside their respective domain S)i and 2)2 to 
establish a common frame of reference in weak isospin 
space. To that end, they need to exchange photons and 
literally see what they get. If they choose to analyze the 
experiment in terms of gauge fields across the boundary, 
they must integrate along the entire path of the photons; 
if they choose to think in terms of the Higgs in the bulk, 
they need only consider the end points of the path. 

We can now answer the question asked at the begin- 
ning of this section: if the VEVs ($) in S)i and 2)2 are 
related by a transformation with U(l) parameter wq and 
SU(2) parameters a pure photon state prepared 
in J)i will be seen by an observer in J)2 as a mix of £{<jj) 
parts photons and 1— ^(a;) parts Z^s and VF^s, with ^(w) 
given by p^ . Physically, photons crossing the bound- 
ary between two electroweak domains are partially con- 
verted to weak vector bosons, which then quickly decay 
to fermions and lower energy photons, leaving us with 
the residual luminosity -^(u;). 

This is the essence of the problem pointed out by Pen- 
rose p5| . If vacuum realignment proceeded at or near 
the speed of light, we would have only 0(1) domain 
boundaries within our Hubble volume, causing glaring 
anisotropics. Unless we are prepared to give up local- 
ity, the resolution of this apparent contradiction between 



^ This can be seen explicitly by diagonalizing the gauge boson mass 
matrix for an arbitrary (<&) |26| 
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standard model and observation lies in the opposite di- 
rection, i.e. in a subluminal rate of realignment allowing 
a larger number of domains to even out such anisotropies. 
(More on this in Section HXj) . 

Since u) was picked randomly at the EWPT, the value 
of can not be predicted. The best we can hope for 
is its probability distribution function (PDF). Viewing ui 
as the SU(2) parametrization 



[no, ni\ cos(aj/2) + ir^— sin(ci;/2) 



of the unit 3-sphere 



1 



(37) 



(38) 



(the Higgs vacuum manifold, up to a factor |(<i>)|) we 
have 



uj/2 = arccos(no) 

UJ Hi 



sin(w/2) 

Substituting Eqs. ([391)- gOl) into Eq. ^ then yields 

£{n) = 1 - 2?!^ sin^ 
i.e. luminosity is determined by 



(39) 
(40) 



(41) 



2 2,2 

n , = + 



(42) 



The luminosity PDF therefore follows from that of 
over the unit 3-sphere. For uniformly distributed n, 
it is simply a step function, Pi{£) > for £ e [1 — 

While the dynamics may modify this result (because 
the electroweak symmetry group is only a subgroup of the 
3-sphere's 0(4)) a uniform n distribution is the null hy- 
pothesis until numerical simulation results become avail- 
able^. 

Having obtained the PDF for one domain boundary 
crossing, Pi{(), we can easily extend our analysis to the 
distribution for the residual luminosity £ = £i ■ £2 ■ ■■■ • in 
after n crossings between n -I- 1 domains with indepen- 
dent ($). By Rohatgi's result for the distribution of the 
product of two stochastic variables [47| . 



P2{i = li-l2)= / dxPi{x)Pi{£/x)/x 



(43) 



and generally 



P„+i(£ ^£i-£2- ... • In) = / dxPi{x)Pn{i/x)/x (44) 




0.05 



0.10 



0.20 



0.50 
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FIG. 6: Log-linear version of Fig. (5] Pn{£) is seen to quickly 
approach a lognormal shape for n > 3. 



TABLE L Mean, median, standard deviation, skewness and 
excess kurtosis of P„{£) (5000 bins). 



n 




£ 


a 


71 


72 


1 


0.7774 


0.7780 


0.1283 


0.0082 


-2.970378 


2 


0.6044 


0.5904 


0.4491 


0.3564 


-2.949819 


3 


0.4700 


0.4550 


0.1362 


0.5845 


-2.944213 


4 


0.3655 


0.3488 


0.1232 


0.7681 


-2.946776 


5 


0.2842 


0.2678 


0.1078 


0.9285 


-2.953191 



(The argument £/x can be shifted around between the 
two PDFs in the integrand by a trivial change of vari- 
ables; for numerics, the form shown is preferable, since 
it minimizes interpolation on the grid at maximum gra- 
dient). R code for Pn{£) is given in Appendix [Cl 

PDFs and cumulative distribution functions for n G 
[1, 5] are displayed in Fig. [5l[n]and[71 Visually, their most 
striking feature is the swiftness by which Pn{i) morphs 
from a step function for n = 1 to a lognormal for n > 3. 
But the real highlights of this analysis are the first two 
records in Table HI they tell us to expect a residual lu- 
minosity £ ~ 0.78 ± 0.13 for photons moving between 
adjacent domains, and £ ^ 0.60 ±0.14 for photons cross- 
ing two domain boundaries, well in line with the £/£e 
ratio of high-z supernovae, ^ 0.79. 



Originally, the PDF was obtained by sampling an lj grid, based 
on the assumption that uii , ui2 and 0^3 are independent, uniformly 
distributed stochastic variables G [0, 4tt] . This skewed the result 
toward higher luminosity, (£) ~ 0.838 vs. 0.778 for sin(0vK) — 
0.22216, but did not affect standard deviation much (cr ~ 0.133 
vs. 0.128). 



VIII. SUPERNOVA DIMMING 

The state of the art of £/£e determination is well il- 
lustrated by 0, where Riess et.al. quote an intrinsic 
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FIG. 7: Cumulative probability distribution functions Pn{L < 
£) for the residual luminosity £ after n boundary crossings 
between n + 1 domains with independent {$) . 

dispersion for individual supernovae (due to sparse sam- 
pling and noisy photometry) of 0.15 magnitudes (in line 
with e.g. ^4,] and [5]), and actual dispersions about best 
fits (on a "gold set" of 157 SNe la carefully selected from 
all available data) of 0.27 magnitudes for 0.1 < z < 1.0, 
increasing to 0.29 magnitudes for z > 1.0. By the stan- 
dard magnitude/luminosity relation 

h/h = lOO^"'''"'^')/^ (45) 

this translates to an intrinsic luminosity uncertainty 
±100°-i5/^ ~ ±0.15, with actual disper sion rang- 

ing all the way from 0.78 to 1.28 for 0.1 < z < 1.0, and 
from 0.77 to 1.31 for z > 1.0. From these figures, mak- 
ing the (intentionally overoptimistic) assumption that all 
errors are independent, the 1 a width for the combined 
e/lE ratio is e w 0.25/7157 w 0.02, i.e. 

f/^B « 0.79 ±0.02 (46) 

In terms of electroweak domains, the statistics is 
poorer. Under the null hypothesis that our ($) domain is 
an average one, it has ^ 14 nearest neighbor boundaries 
(see Section |V|. With high-z supernova searches concen- 
trated to the equatorial plane (so that the evolution of 
light curves, needed for distance determination, may be 
tracked from observatories on both hemispheres) we are 
probably looking through fewer than half of them (see 
the section drawing in Fig. [3]). Based on Table U the 1 a 
expectation width for one boundary crossing is therefore 
roughly 0.128/%/7 « 0.05, i.e. 

~ 0.78 ± 0.05 (47) 



This suggests that there is a single layer of domain 
boundaries between us and high-z supernovae^. 

The good match with one domain boundary allows a 
consistency check: the boundary distance curves of Fig. 
|4] illustrate the fact that dimming should not start too 
close to home, lest we end up expecting too many domain 
boundaries and too dark a universe at higher redshifts. 
By inspection, the "safe" redshift is z > 0.3. Ideally 
(assuming that we are near the center of our domain) we 
would like to see no sign of dimming at all for lower z 
than this. 

This condition is satisfied: in their 2003 analysis of 
published SNe la data, Padmanabhan and Choudhury 
[4^ found that it split neatly in two subsets about z = 
0.25, consistent with different luminosity zero points and 
no continuous dimming at all - exactly what crossing 
of a single domain boundary at z = 0.25 should look 
like. Their foUowup analysis with even more data (49j 
showed that the split point can be moved all the way to 
z = 0.34 with little consequence, definitely bringing it 
into our "safe" zone. The conclusion that there is little 
or no evidence of accelerated expansion for z < 0.3 was 
again reached by Shapiro and Turner in [H^l ■ 

Summing up, the SNe la data appear to fit our expec- 
tations well, based only on the standard electroweak and 
big bang models and without any need for a cosmological 
constant. 



IX. DISCUSSION 

This paper presents a novel picture of the universe as 
a foam-like structure of large-scale electroweak domains 
in a scaling state, arguing that this is both a natural 
consequence of the standard big bang and electroweak 
models and a good match to SNe la observations - in 
fact the only one to date sans A > 

In this picture, "dark energy" (in the restricted sense 
of an unseen contribution to Q) is carried not by a hy- 
pothetical new (set of) scalar field(s) with negative pres- 



The higher (£) of the original PDF left some room for multi- 
ple boundary crossings, with a best fit including one- and two- 
boundary crossings a 65/35 mix. Assuming that an effective 
cutoff distance rcut can be defined, such that observations are 
made in a sphere of radius rcut, this suggested an average dis- 
tance to the second domain boundary ^ 0.65^^/^ rent ~ O.Qvcut 
and led to a very rough guesstimate of 2 ~ 1. 
After completing the first version of this paper, the author be- 
came aware of a dimming model which had been proposed shortly 
before that version was written _54|. In spite of attractive phe- 
nomenological similarities - a domain structure causing dimming 
in discrete steps - the underlying physical motivation is quite 
difl'erent: 54] invokes a stringy braneworld scenario with com- 
pact extra dimensions to create domain walls where photons mix 
with a hypothetical "para-photon" , leaving unresolved the prob- 
lem of reconciling the required low energy density of the walls 
with experimental limits on the string energy scale, already in 
the hundreds of GeV. 
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sure, as in the popular quintessence models, but rather 
by the long wavelength modes of a very familiar one: the 
standard model Higgs (or the superset of Higgs fields em- 
ployed by standard model extensions) . The SNe la lumi- 
nosity deficit is caused not by accelerated expansion, but 
by photon conversion at domain boundaries. No exotic 
new physics is needed. The most "speculative" assump- 
tion is that electroweak domain growth is an unexcep- 
tional case of generic 3D domain growth. 

In the reader's mind, this picture has likely already 
been met with a strong objection (if not earlier, then 
upon review of Fig. [4]): what about the CMB? If we 
live inside an irregular polyhedron with faces of differ- 
ent opacity, should we not see huge anisotropics in the 
microwave background? 

The first part of the answer is a counterquestion: 
what's a poor photon to turn into? The 2.73 K of the 
present-day CMB translate to 2.35 • 10~'^eV, well be- 
low the rest mass of any known particle except the pho- 
ton and - perhaps - the lightest neutrino. With cur- 
rent neutrino oscillation data indicating mass differences 
> 10~^ eV [511 (plus a much debated claim from the 
Heidelberg-Moscow double beta decay experiment for 
nil, > 0.17 eV (Hli) conversion of CMB photons may sim- 
ply be kinematically forbidden, making electroweak do- 
main boundaries effectively transparent to them. In this 
regime, instead of conversion, photons can be expected to 
undergo random scattering, resulting in the diffuse glow 
seen in ordinary foams. 

Two bounds on the lightest neutrino mass are thus 
implied: rrii, > 10~^ eV to prevent readily visible CMB 
anisotropies, rrij/ < 1 eV for supernova dimming with- 
out the reddening problem of dust models (the optical 
range being roughly 2 to 3 eV). Incidentally, the lower 
bound also explains an observed discrepancy in the rela- 
tion between luminosity and angular diameter distances 
for supernovae and radio sources which has been used 
to argue against photon loss as an alternative to Hubble 
acceleration [5^ : if domain boundaries are transparent 
to microwaves, they are necessarily transparent to radio 
waves, too. 

It may be possible to improve the lower bound on 
by going back in time, to higher CMB temperatures. If 
'^1' ^ 0.13 eV, Heidelberg-Moscow notwithstanding, we 
will eventually reach a point along the road to recom- 
bination at 0.26 eV ^ z ^ 1089 where CMB photons 
could be converted to neutrinos (the ionization poten- 
tial of hydrogen is 13.6 eV, but recombination occurs 
later due to photons in the blackbody distribution tail 
and transitions between excited states). Two consider- 
ations then come into play. First, the number of do- 
main boundaries cutting through the equatorial plane at 
a given distance is roughly proportional to the number 
of domains out to that distance. The WMAP angular 
resolution, « 0.3°, is enough to resolve 1200 individ- 
ual faces, equivalent to ^ 1200/6 = 200 domains. Any 
N,j, > 200 • 0.26/(2m^) should therefore be safe from di- 
rect observation by WMAP. Given rrii, < 0.13 eT^, this 



implies > 200. It then becomes a question of statis- 
tics: how large are the anisotropies caused by differences 
between (individually unresolvabe) faces? 

Fluctuations in primary radiation which was exponen- 
tially damped by O(N^) boundary crossings are clearly 
irrelevant: > 200 implies a residual luminosity {£) < 
0.782°" - 10-22; with = 10^ as suggested by Section 



lYl (i) 
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Rather, the dominating contributions 



should be from secondary photons which were redshifted 
below the 2m,y threshold after only a few boundary cross- 
ings. Now, Fig. [6] reminds us that residual luminosity 
distributions become approximately lognormal, implying 
constant relative standard deviation, already at n = 3. 
To estimate the resulting anisotropies, we can therefore 
simply read off cr/ {£) ~ 10"^ from TableHl invoke the cen- 
tral limit theorem and divide by the square root of the 
number of faces covered by a pixel at the 2m,y threshold, 
A^o ~ (2toi,/0.26) • Ng, domains out. 

For a sphere with radius A^Oi the area of a cap sub- 
tended by angle 9 is 2nN§{l - cos(6i/2)). Tiling it with 
tetrakaidecahedra of unit diameter gives us 7/(7r/4) faces 
per unit area, for a total of 28A'o (1— cos(0/2)) faces. This 
is actually an underestimate, since it does not take the 
uneven spacing of domain boundaries into account; at 
large z, we should rescale the radius by the comoving 
coordinate distance 



D{z) 



2Q.Z +{VL- 2)^1 + Viz - 1) 



n'^{i + z) 



(see e.g. [55|) to get 



N K2d, {NQD{z)f (1 - cos(6'/2)) 



(48) 



(49) 



(another way to see this is to first correct tile size for 
angular diameter distance and then Aq for scale factor). 
Inserting the WMAP angular resolution 9 — 0.3 • 7r/180, 
n = l, z = 1089 • 2mi./0.26, my ^ 0.1 eV and = 10^ 
yields fluctuations ~ 1Q~^/^/N ~ 10"^, in line with ob- 
servation. For smaller values, reduce Q. (everything else 
being equal, the old preferred SI « 0.3 buys a further 
factor 1/10) and/or increase rriy and/or (by equation 
(I32p there is plenty of margin here too; the order of mag- 
nitude argument pointing to A^$ ^ 10^ is certainly very 
rough) . 

While this estimate is already in the right ballpark, 
keep in mind that it ignores smoothing by random scat- 
tering (diffuse glow) and the successive injection of high 
energy radiation from astrophysical sources. Over time, 
down-conversion and scattering of this radiation by do- 
main boundaries should also build up a contribution to 
the CMB. Intuitively, its spectrum should be the unit 
eigenstate of a Markov chain through the decay channels 
available at each photon's energy, with a peak about the 
average energy of photons produced by the last decay in 
the chain; essentially the "resonance peak" of the light- 
est particle capable of decaying to photons (yes, there is 
a definite iconoclastic possibility lurking in this observa- 
tion). Computing it should not be hard, given a reliable 
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Monte Carlo generator of highly virtual weak boson de- 
cays (alas, apparently not a priority in off-the-shelf simu- 
lation packages geared toward high energy experiments) . 

A detailed CMB calculation would also require all neu- 
trino masses to be specified; the a and (3 parameters of 
Section lVTl to be extracted by numerical simulation of the 
scaling state in the electroweak model under considera- 
tion, in turn requiring knowledge of the Higgs mass(es); 
and of course a consistent choice of cosmological expan- 
sion history. It would be a large undertaking, but a 
systematic exploration of model-|-parameter space might 
yield useful constraints on candidate cosmological and/or 
extended electroweak models. 

This resolution of the Penrose conundrum 45] has a 
problematic implication: the CMB is not quite the fossil 
we thought it was. Long after it left the surface of "last" 
scattering, it was interacting with domain boundaries, 
and later on it may have picked up significant (dominat- 
ing?) contributions from other sources. This begs the 



question if it really is telling us so much about the early 
universe as we currently like to believe. Its large scale 
isotropy would undeniably be difficult to explain differ- 
ently, but at smaller angular separations, its features may 
be of more recent origin. In particular, while none of this 
directly contradicts inflation (even 51 « 1 seems possible 
to accomodate, at least at this early stage), it does call in 
question the relevance of the CMB towards testing both 
inflationary and alternative scenarios [H^ [s^l ■ 

On the other hand, once the large-scale electroweak 
foam picture starts sinking in, many new venues of inves- 
tigation readily suggest themselves. Do domain bound- 
aries balance the energy budget of the universe? What 
are their effects on structure formation? Can the evi- 
dence for dark matter from cosmic shear be reinterpreted 
in terms of domain boundary effects? Can they help ex- 
plain observed CMB anomalies ^58] ^59] [60.]? 

The answers are out there. 
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znVector <- array(NaN, c(l, maxN)) 
. php? a5t'l'8ig^ip'ciftfgg2a&id=81 ) . 

J. ± uiiu L±uii V, L y \ 

alpha * sqrtd + beta2*(t - t_alpha)) } 

integrand <- function(t) { 

0.5 * alpha * beta2 * a(t) / 
sqrtd + beta2*(t - t_alpha)) } 

rprop <- function(t) 

{ 

integral <- integratedntegrand, 
t_alpha, t, 
rel.tol=lE-13, 
abs.tol=lE-13) 

(a(t) * r(t) + integral$value ) / a_alpha 

} 

dn <- functionCt, n) { (2*n - 1) * rprop(t) } 

# Derived parameters 
a_alpha_dev <- f unction (t_alpha) { 

(a_alpha - a(t_alpha) ) "2 } 
t_alpha <- optimize (a_alpha_dev, 

c(0, 1), tol=lE-30)$niininium 

alpha <- hR(t_alpha) * r_alpha 
beta <- 2 * v_alpha / alpha 
beta2 <- beta*beta 

# Main loop 

t_em_dev <- functionCt, n) { 
(dpropCt, 1) - dn(t, n))"2 } 

forCn in l:maxN) 
{ 

t_ein <- optiinize(t_ein_dev, 

c(0, 1) , n = n, tol=lE-30) 
z <- a(l) /a(t_eni$minimum) - 1 
znVector [n] <- z 

if ( (t_ein$ininiinum < t_alpha) I I 
(z > maxZ) ) break 



znVector 



APPENDIX A: DOMAIN BOUNDARY 
POSITIONS (R CODE) 



# Main program 
redshif ts () 



# Replace flat dust with pet cosmos here 
w <- 

q <- 2 / (3 * (1 + w)) 
a <- function(t) { t"q } 

# Auxiliary cosmological functions 
arecip <- function(t) { l/a(t) } 
dprop <- function(tO, tl) { 

integrate (arecip, tO, tl, 
rel.tol=lE-13, abs.tol=lE-13)$value } 
hR <- function(t) { a(t) * dprop(0, t) } 

# Compute expected domain redshifts 

redshifts <- functionC 

a_alpha = 1/1090, # a(t_alpha) 

r_alpha =1, # r(t_alpha)/rH(t_alpha) 

v_alpha = 0.00022, # v(t_alpha) 

maxZ = 5, maxN = 10) 



APPENDIX B: LUMINOSITY DISTRIBUTION 
FOR N = 1 (R CODE) 

sw2 <- 0.22216 # sin"2(theta_W) 

oneCrossLuminosity <- function(n) 

# Return PDF and CDF of luminosity 

# Simple step function for input to nCrossLuminosity () 
{ 

pdf <- array (0, c(n)) 
cdf <- array (0, c(n)) 

i <- round(n*(1.0 - 2.0*sw2)) 

pdf[i + l:(n - i)] <- 1.0/(n - i + 1) ; 

cdf[l] <- pdf [1] 
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for (1 in 2:n) { cdf [1] <- pdf [i] + cdf [i-1] } 

# Return results 

list (pdf =pdf, cdf=cdf) 

} 

# Main program 

sample <- sampleLuminosity (5000) 
save(sample, file=" sample . txt" , ascii=TRUE) 



TABLE XL Parameters of fit and error squared for CDF 
parametrization F{£;a,b,c,d) (|Df |) . 



n 


a 


b 


c 


d 




i 


0.29650820 


0.4238114 


-0.1885901 


0.05361303 


6.011231 


2 


0.05948206 


0.7159291 


-0.3665269 


0.1868697 


-0.03652571 


3 


0.09153357 


-2.05914 


72.6966 


14.41858 


-0.3107125 


4 


0.01036675 


0.3067824 


-0.3199721 


0.1081582 


-0.08748841 


5 


0.4831496 


-2.23015 


31.56453 


22.18672 


-0.05052748 



APPENDIX C: LUMINOSITY DISTRIBUTIONS 
FOR N > 1 (R CODE) 



nCrossLuminosity <- function (maxCross, nlpdf , nlcdf) 
# Compute PDF&CDF for 2 to maxCross boundary crossings 
{ 

n <- length (nlpdf ) 

pdf <- array (0, dini=c (maxCross , n)) 
cdf <- array (0, dim=c (maxCross , n)) 

for (1 in l:n) 
{ 

pdf [1, i] <- nlpdf [[i]] 
cdf [1, i] <- nlcdf [[i]] 

} 

if (maxCross > 1) 
{ 

for (nCross in 1: (maxCross - 1)) 
{ 

cumulative <- 
for (iy in l:n) 

{ 

integral <- pdffnCross, n] * nlpdf [iy] / iy 

if (iy < n) 
{ 

for (ix in (iy+1) :n) 
{ 

integral <- 
integral + 

pdffnCross, (n*iy) '/,/'/, ix] * 
nlpdf [ix] /ix 

} 

} 

integral <- n*integral 
pdf[nCross + 1, iy] <- integral 
cumulative <- cumulative + integral 
cdf[nCross + 1, iy] <- cumulative 



} 



} 



# Return PDFs and CDFs 
list (pdf =pdf, cdf=cdf) 



APPENDIX D: DISTRIBUTION FITS 

Forn > 1, the CDFs P„(i < £) of Fig. [7] are found to 
be well parameterized by 

F{i; a, b, c, d) = , ^ ryr- (Dl) 

F,„(l;6,c) exp(d) 

where F\n{x; b, c) is the lognormal CDF 

F,„(.;^c).i{l + erf(i^^)} (D2) 



erf(x) is the error function 



2 /"^ 

erf(a;) = —j= I dy exp(-?/^) 



(D3) 



and the parameters a,b,c,d (plus total error squared e^) 
obtained on a 5000x5000 point grid are given in Table 

IE 

Differentiating (|Dip in £ yields the corresponding 
parametrization of P„ (£) : 



/(£; a, b, c, d) = a, b, c, d) 



b, c) + b, c)) exp(dr ) ax 

Fin{l;b,c) exp(d) 

where fin{x; b, c) is the lognormal PDF, 



(D4) 



1 / (lna;-6)-'\ , , 
fin{x; b, c) = = exp (D5) 



} 



2cxV27r ^ \^ 2c2 y 
For the first couple of n, fitting directly on P„(€) with 

g{e; a, b, c, d, e) ee a, b) | Y~~j: + ^^ + 62:1 (D6) 

achieves a better reproduction of the sharp peak domi- 
nating the PDF. Parameters obtained on a 5000x5000 
point grid are given in Table Hill 



# Main program 

multi <- nCrossLuminosity (5 , sample$pdf, sample$cdf) 
save(multi, f ile="multi .txt" , ascii=TRUE) 
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TABLE III: Parameters of fit and error squared for PDF parametrization g{£;a,b,c,d,e) (ID6|) . 



n 


a 


L 



c 


a 


e 


J2 

e 


1 


2.4843356-05 


-0.1780371 


0.2699788 


5.571233e-08 


0.0002125698 


8.1680096-05 


2 


1.1985346-06 


-0.5352473 


0.2612992 


-1.576451e-07 


0.0001519190 


9.1032216-05 


3 


1.8979276-07 


-0.6941752 


0.332284 


-2.835761e-08 


0.0003936519 


-0.0003655848 


4 


8.4990976-08 


-0.9647373 


0.3736279 


-1.269617e-08 


0.0003422317 


-0.0003472561 


5 


5.9192396-08 


-1.235350 


0.4105253 


-1.767185e-08 


0.0003091985 


-0.0003447646 
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